Non-planar double-box, massive and massless pentabox Feynman integrals in negative 

dimensional approach 

Alfredo T. Suzuki*, Alexandre G. M. Schmidt 

Universidade Estadual Paulista — Instituto de Fisica Teorica R.Pamplona, 145 Sao Paulo - SP CEP 01405-900, Brazil 

(February 1, 2008) 



O 
O 

(N 

> 

o 
O 



X 



Negative dimensional integration metliod (NDIM) is a technique which can be applied, with 
success, in usual covariant gauge calculations. We consider three two-loop diagrams: the scalar 
massless non-planar double-box with six propagators and the scalar pentabox in two cases, where 
six virtual particles have the same mass and in the case where all of them are massless. Our results 
are given in terms hypergeometric functions of Mandelstam variables and for arbitrary exponents 
of propagators and dimension D as well. 
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I. INTRODUCTION 



Studies in particle phenomenology require more and more sophisticated calculations and the measurement of 
I ■ the (g — 2) factor has now an error of 1 ppb order of magnitude thanks to the perturbative approach. Within 
^ this perspective there have been some interest in massless double-box — planar and non-planar |^ — and pentabox 
[ — integrals. Smirnov jj] studied the former, using the Mellin-Barnes (MB) technique, in three cases: scalar and tensorial 
[ with four legs on-shell and scalar with one leg off-shell; Tausk and Smirnov and Veretin [Q], also using the 
■ MB method, presented explicit results for non-planar, or crossed, double-box with seven and six propagators, in 
' the special case where the exponents of propagators are equal to one; Kiimmer et al came across with the same 
integrals studying the potential between quarks in the Coulomb gauge |Q; Anastasiou et al calculated the latter, in 
the integration- by-parts approach, for both the scalar and tensorial cases with massless internal particles Q. The 
results which are obtained with the MB approach are, like the ones obtained by NDIM, expressed in terms of infinite 
• series of hypergeometric type. Of course, progress along this line is greater in covariant gauges, however perturbative 
I ' calculations in non-covariant ones are also carried out and need sometimes more powerful techniques than the formers 

(-H I In this paper we choose to tackle the scalar on-shell non-planar double-box integral for arbitrary exponents of 
• • . propagators, a result that is missing in the literature and the pentabox integral in two cases: where all particles 
are massless, and where six virtual particles have the same mass, a diagram which, for instance, contributes to 
photon-photon scattering, and as far as we know was not yet calculated before. NDIM gives us several results in 
\^ terms of Mandelstam variables and masses, each valid in a certain kinematical region. We give results in terms of 
9^ ' hypergeometric series for arbitrary exponents of propagators and dimension. In our approach no reduction formulas 
or integration-by-parts methods are used or even necessary. It is also worth observing that NDIM is a technique 
which can be applied to other gauge choices like the Coulomb and the light-cone gauges ||ll|. 

An important feature of NDIM is that it is not a regularization technique. It is worth remembering Dunne and 
Halliday ]l^ : the negative-dimensional integrals (in Z3-dimensions) can be related to positive dimensional ones (in 
2iV-dimensions) over Grassmannian variables; in fact, one has just to make D ^ —2N. So, in NDIM context there 
are no singularities, no poles etc. However, when we perform the analytic continuation, in order to allow negative 
exponents of propagators and positive dimension, then poles appear for specific values of those exponents and physical 
D — A dimension and we have the same results which other techniques provide. This is therefore a consistent method 
to solve Feynman loop integrals pertaining to usual covariant or to non-covariant algebraic gauges, like Coulomb and 
light-cone ones (even at the two- loop level). 

The aim of our paper is not to establish the axiomatic foundations for NDIM nor to demonstrate in a rigorous 
mathematical sense its principles and basis. Instead, given the simple steps that the methodology requires to work 
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out complicated Feynman integrals we are interested in testing it to the limits of our present calculational abilities. 
For this purpose we are presenting here in another true two-loop calculation the exact results yielded by this 
method. Such results must be compared with others obtained using different techniques so that not only previous 
answers be double-checked and the confidence in the novel method be increased but also in order to demonstrate the 
feasibility of the latter in performing the calculations. 

The outline for our paper is as follows: in the next section we study pentabox integrals, where in one case virtual 
particles are massless and in the other case, we consider six of the virtual particles massive with equal masses, a graph 
which contributes to photon- photon scattering (which were considered recently in [p^ ), although a full calculation 
of such an effect as well as evaluation of beta functions in physical processes are beyond the scope of the present 
work. We solve also the scalar massless non-planar double-box with six propagators. In section 3, we present our 
conclusions. 

II. PENTABOX AND NON-PLANAR DOUBLE-BOX INTEGRALS 

Let us define immediately the relevant three negative-dimensional integrals, namely, 

= 1 1 d''qd%-P{q,h,p,p',Pi), (1) 

^^"/ /^^^c?^ MP(9,fci,p,p',pi,/i), (2) 

and 

MP6 = J J d''qd''rNPeiq,r,p2,P3,Pi) , (3) 
where the integrands are respectively 

P ^ {qy{q-p)'Hq - - fcl -Plfiklriq - Plf'^iq-p- pT : (4) 

MP ^ (g2 - ^2). [(^ _p)2 _ ^2]. [(^ _ ^^)2 _ ^^2]'c _k,^p,f~ ^i'Ukl^ (5) 

X [{q-p,f-^^X[i1-p-p'f-^^"Y ^ 

NPe ^ (g')'(g -P3)'^('Z + rf\q + r+p2Y\r^r{r-piY- , (6) 

which represent massless pentabox, massive pentabox and non-planar double-box, respectively. Once they are intro- 
duced, we evaluate them using NDIM. 

When Halliday and co-workers advanced the idea of negative-dimensional integration they proved that it is equiv- 
alent ||l^ to Grassmannian integration in positive dimension, the correspondence being as simple as D ^ —2N. This 
fact is implied in the very structure of the above integrands. 
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A. Massless Pentabox integral 



q — r 



pi 
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FIG. 1. Scalar massless pentabox with all external legs on-shell. Mandelstam variables are defined as s = {p + p')^ and 
t=ip~pif 

The negative-dimensional integral for massless pentabox diagram — see Figure 1 — has the following generating 
functional, 



2\ -D/2 



-- {aX2UJS + Pr]X2t + f]0(t>t) 



(8) 



where A2 = 7 + 6* + and A = {a + /3 + rj + uj)X2 + 7(/) + 0(j). 

Taylor expanding (W) we have our negative-dimensional integral "P as a factor in a seven-fold summation series, 



i.j,k,l,m,n,r—0 



i\j\kU\m\n\r\ 



(9) 



On the other hand taking (H) and making an expansion (including a multinomial one) in power series, we obtain, 



^123 +^123 ^^4567+^456 rv^l4+^l 47 + -2'l3 /9^257 +^^258 + -2'24 ^^367 + ^369 + 2'i 25 7^^'456 + Vygg ^ ^ jXi23 + -2'34, 



?7' 



Xi\...X7\Yi\...Yq\Zi\...Z5\ 

X (-s)^i^^ (-X1234567 - D/2% 



X,Y,Z=0 



(10) 



where sum over X,Y,Z is a shorthand notation for a 7-, 9-, and 5-fold sums respectively. Moreover, we define 
X12 = Xi -I- X2, so on and so forth. 

Therefore, the exponential above generates a series indexed by 21 indices, while the 7 propagators in the argument 
of the integrand give rise to 7 equations and the multinomial expansion to another one, see eq.(j9|). Now, solving 
both equations for V we conclude that there must be some relations among the two sets of indices {X, Y, Z} and 
{i,j, fc, I, TO, n, r}. It is a system of algebraic equations. 
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for which there is not a unique way to solving it, since we have 21 "unknowns" and just 8 equations. In fact, there are 
203,490 possible 8x8 systems which can be solved in terms of exponents of propagators z, j, k,l,m, n, r, dimension D 
and some of X, Y, Z. The computer can calculate such 8x8 determinants easily: 134,890 of them are zero, i.e, give 
empty sets of solutions. 

So, our result will be written in terms of a 13-fold series of hypergeometric type. This can be worked out conveniently 
and be simplified, since we have three possible variables: t/s, s/t, and unity. Series which depend on both t/s and 
s/t can not be convergent, so we will not consider them. 

Next, our strategy is to search for the simplest hypergeometric series among the remaining 68, 600 solutions. The 
criterion for this search is dictated by the fact that the more sums with unity argument the simpler it is, since one can 
sum them (at least in principle, respected some relations amongst the parameters) using Gauss' summation formula 

2Fi(a,6;c|l) = ^M£_^, Re(c - a - 6) > 0. (12) 
1 (c — a)l (c — 0) 

In other words, when some of those 13 sums can be rewritten as gamma functions. Of course, not all sums with 
unity argument can be summed this way, as we will shortly see. These cases occur when the resulting function is of 
the type 3_F'2(---|1), or even more complex ones. 

There are many different ways in which the 13-fold series appear, and we can classify them according to the following 
general form. 



(Momenta) (Gammas' 



00 ^ 

y f 

Til!... rig Itoi!... 7714! 

all— 



where z is any one of the three possible variables t/s, s/t, or 1, and A represents either ni23 = ni + n2 + 713, or 

«1234 = J2i=l "ii or ni2345 = J2i=l '^i- 

The simplest hypergeometric series representations for the scalar massless pentabox is given by hypergeometric 
series with three such "variables" (here we use the word "variables" , to denote the one variable s /t appearing thrice 
as summation variables within the series): 

= s^-^ t' Pi"" E(Y4, Y5, Ye) 3i^2(a3, ^3, C3; 63, Ml) (13) 

where 

{-k\Yi){] -k-l-m-n- D/2\ - Yj.eKD + k + l + m + n + rlY^^e) 
2.iY4,Y5,Y6j- 2^ ^ k - m - D /2\ - Y,)iD + k + I + m\Y,,,) ^'^^ 

{D + i + k + l + m + n|r456)(fc + I + D/2\Ye){-j\Y^56) (sA)^«'= 



(l + a~j|y456) YilY^lYel' 



and the parameters of 3^2 are 



a3 = -fc + y4, (15) 
63 = j — k — I — jn — n — D/2~ Y^^g , 
C3 = -k - I - m - D/2, 
63 = 1 - fc - m - D/2 - Y5, 
f3 = -k-l-m-n- D/2 



with a = i + j + k + l + m + n + r + D and 
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P^'^ = {-i\a~j){~k-l-m-n-D/2\j){-l\k + l + m + D/2){-m\-k-l-D/2) 

X {-r\k + l + r + D/2){(T + D/2\ ~ 2(j - D /2 + j)(k + I + m + D\i + n) (16) 

In the above equations, we have used the subscript "3" and superscript "AC" to mean that we have three "variables" 
and the result is analitically continued into positive dimension D. Pochhammer symbols and some of their properties 
are used throughout the expressions, 

(a|6) ^ (a), ^ ia\-k) = jt}l_^ {a\b + c) = (a + b\c){a\b). (17) 

Observe that in (|l^) there is a fourth sum, namely in Zi, which has unity argument and could not be summed 
because it is a 3F2 hypergeometric series which can be expressed in terms of gamma functions only in some special 
cases. This means that we were able to sum up up to 9 series (with unity arguments) using (^|). 

The second type of result provided by our method is given by hypergeometric series with four "variables" , 



i-j\Xi567){-n\X456){-k\Xi)i-l\X57)i-Cj\X4567) 



^^^^ [i - CT|X4567)(r - Cr|X4567)(-fc - l\Xi57) 

^ {D/2 + m\Xi5){D/2 + k + l\XQ){-k ~ I - m - D/2\Xr) (Vg)^^^°^ 

^ {D + k + I + m\X456) X4lX5lXelXjr ^ ' 

where the subscript "4" means a 4- fold hypergeometric series representations for V, and 

P/'^ = {-i\cT){-r\a){-k -l\-m- D/2){a + D/2\ -2a- D/2){k + l + Tn + D\-m- D/2){-m\2m + D/2). (19) 

Note that this four-fold series can be used to study forward scattering. Taking i = we arc left with only the first 
term in the series, which is equal to unity, that is, the series collapses and the result is merely (the superscript FS for 
forward scattering case) 



The last one, given by hypergeometric series with five "variables" , 

-p^ _ ^i+j+r+D/2 ^k+l+m+n+D/2 pAC 

{-n\Y^^c,){-k\YJ + Zi)(-fc -l-m- D/2\Z^2){k + I + D/2\Y^) 

Y z=o i'^-k-m-D/2\Yrs9){l+J-k-l-m-n-D/2\Y7sg + Zi2){D + k + l + m\Y7S9) 



^ ij_±r + D/2\Y7S9 + Zi2){i+J+D/2\Yjs9 + Z,2) {s/tY^-^> 



(20) 



where 



{l + i+j+r + D/2\Yn9 + Zi2) Y7\Ys,\Y9\Zi\Z2\ ' 

Fg4C ^ (-_-| -j-r- D/2){-j\k + l + m + n + D/2){k + l + Tn + D\ - 7n-D/2){-l\ -k-m- D/2) 

X {-m\i + j + 'm + D/2){-r\k + 'm + r + D/2){a + D/2\j + r-a). (21) 

Observe that in the last line of ( pO| ) when we take either i = — lorr = — 1 those pertinent Pochhammer symbols 
within the series do cancel out, simplifying it. 

Of course, there are results which depend on the inverse of such variables, i.e., V^it/ s), VAis/t) and V5{t/s), which 
means an interchange of one or more pairs of exponents of propagators and s <-> 

How are these different expressions related to one another? How do they relate to previous calculations by other 
methods, for example for special cases of the propagator powers? How unique or ambiguous is the analytic continuation 
in dimension to get back to answers in positive dimensions? What form of renormalization is necessary to connect 
with conventional calculations? 

To answer the first question one must recall that if two series {Va and , for instance) represent the same function 
(the integral V) then they must be related by analytic continuation |T^. So, the results provided by NDIM are always 
related by analytic continuation either directly (overlapping regions of convergence) or indirectly (no overlapping of 
regions). In the present case there are no previous calculations in the literature considering arbitrary exponents of 
propagators. 
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B. Massive Pentabox integral 



Introducing masses in the NDIM context is very simple |T^. The generating functional becomes, 



G 



MP 



d^'q d^'ki exp {-a{q^ - Ai') - m - pf - /i'] - 7[(g - k^f - /?] - e[{q ^ k, ~ pif 
-c^kl - Ti[{q -pif- -u:[{q-p~ p'f - ^?] ] 



2\D/2 



exp 



A 



(Q;A2a;s + I3r]\2t + P9(f>t) 



exp[{a + P + j + e + ij + uj) 



(22) 



(23) 



so that we have six additional sums, which are generated by the second exponential in the second line above, which 
corresponds to the massive sector. Besides, the resulting hypergeometric series will have variables of the form. 



t 



their inverses or unity. Powers of them also occur, like \fsjt and {s/tY. 

Once more we look for convergent series. Among the very BIG number of possible systems — altogether 
27!/(8!19!) ~ 2,220,075 — there are 1,093,289 which have no solution, and the remaining 1,126,786 among which 
we are able to find solutions. So we are left with 49.24% of the total, from which we hunt for the most convenient 
ones. 




FIG. 2. Pentabox diagram with six massive propagators. We consider the case of equal masses and external particles on-shell. 
At two-loop level it contributes to photon-photon scattering. 

The simplest hypergeometric series representation is given by a seven-fold summation, 

(MV\ = ^( pAC (-»l^l23)(-j|^4567)(-fc|A:i4)(-?|A:257)(-'-»|^367)(-"|^456)(-r-|A:i23)(-a|Xl234567) 

^ " ' ' (-rn-<T|Xi245 + 2X367)(-fc-Z|Xi2457)(i)/2|Xi234567) ^ ' 

{-i-j-m-n- r - D/2\Xi2i57 + 2A:36)(-/c - I - m - D/2\X-j){Dl2 + m|Xi245) ( s \ ( t ^ ^^'^ 

X 



(-Z- j-n-r|2Xi23456 +^7)^i!-^7! VmV W 
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where 





P,''^ = {D/2\m){-a\-m){-' 


i — j—n — r\- 


- m 


-D/2){-k-l\ - 


- m — 


D/2). 




(25) 


As a sample numerical calculation we give in 


the Table 1 below 


an expansion in 


the e : 


= 2- 


D/2 parameter. 


TABLE I. 

wcYC suimii('( 


Column "terms" indicates whore we truncated all the 
1 we get a good precisitni.Tlie result is (jf tlit^ forui .1 -i 


series. Observe that even for 
-Bf-Cf- where A.B.C are 


a few terms in the 
giA't'ii ill tli(> table. 


series that 


Terms 








e 








'2 

e 





0.167027110 






-0.012614420 








0.3914896 


1 


0.167553265 






-0.012917188 








0.3927286 


2 


0.167558252 






-0.012920469 








0.3927776 


3 


0.167558333 






-0.012920522 








0.3927779 


4 


0.167558335 






-0.012920523 








0.3927778 
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Hypergeometric series with ten summation indices, nine "variables" also occur, 

[l-k-l-m- DIW34 - Xi23)(fc + m + i:)/2|X2 - 1^3 - ^i)(-fc -l-m- D/2\Wzi + Zi) 

(1 - fc - ; - Dl2\W^i ~ X3) 

(i + J + r + i?/2| - - X123 + Zi)(l - d - i?/2 1 14^123456) (sA)^^^^ (/iVi) ^1=3456 



(26) 



{-k-l-m-n~ D/2\W3i5 + Zi) Xi\...X3\Wi\...Wq\Zi\ 



where 



Pg^^ = {-k-l-m-n- D/2\n){a + £1/21 -2a- D/2){k + l + m + D\-m- D/2) 

X {-j\a){-l\k + l + m + D/2){-7n\i + j + m + r + D/2). (27) 

Observe that there is a tenth series with unity argument, which is not summable since it is a 3^2(0, b, c; e, /|1). 
Finally the last result we present for the massive pentabox is a 11-fold series, 

{MV)ll = ^^2y+j+r+D/2 ^k+l+m+n+D/2 pAC 

-*|Xi23)(-A:|Xi + W3 + yr + ^i)(-«|W^5 + >789)(-r|Xi23)(A: + I + L'/2|X3 + yg - VK34) 



X 



E 



^jj^Q (1 + j - fc - ; - m - n - i?/2|M^345 + ^789 + ^12) 



{\-k-l-m-n-D\- X123 + ^789 + Z12 + Vt^345)(l - fc - / - m - i^|Ty34 - ^123 - ^789) 
^ (l-k-l-m-n- D\W345 - ^123)(1 - A: - m - D/2\W3 - X2 - Ys + Zi) 

i-i -j-r- D/2\Xi23 - ^789 - Zi2){-k ~l-m- D/2\W3i + Zu) 



X 



{-i-j + k + l + m + n-r + D/2\ - W345 + 2X123 - ^789 - Z12) 



(28) 



where 



Xl\X2\X3\Y^\Yg}.YQ\Zl\Z2\W3\W4,\W^\ ' 

P^^^ ^(^-i-j + k + l + m + n-r + D/2\i + r){k + l + m + D\-m~ D/2) 

X -i-r- D/2){-l\ -k-m- D/2){-m\k + 2m + D/2) . (29) 

It is important to note that all the series are valid within their regions of convergence [|l5|, e.g., in the first one ( [2^ ) 
one must have |s//i^| < 1 and < 1; ''^^'^ so on for others. 

We mention that there are also other hypergeometric series, for instance, 

aU=0 ^ \ ^ all=0 ^ \ / 

and also 10-fold series, with 8 "variables", which has poles in the exponents of propagators, like T{i — r), 

(MP)>-r(i-rn^.)j:[ ) +r(,-,)r(,Ji: iL) , ,31, 

all=0 / \ / an=0 ^ \ / 



which is regularized introducing |18yi9f| for example, i = —1 + S and then expanding around (5 = 0. 

Among all the possible (more than a million) hypergeometric series representations for the integral in question there 
must be many relations among them by analytic continuation (directly or indirectly). Yet, the study of them all is 
such a formidable task that for practical purposes it is virtually impossible to even think of doing it thoroughly, so 
we can only conjecture that the entire set of hypergeometric series representations cover all the kinematical manifold. 
Also, these analytic continuation formulas can in principle be obtained from them. 

It is worth remembering some important points in the complex analysis theory |l^ ]. The reader can observe that 
some of the results have poles of different order than others, i.e., some of them are like l/e" and other(s) like 1/e^, 
with a ^ b. In our previous work on box integrals for photon-photon scattering the same occured ]l9[ |. In that case 
we had a branch cut at s = 4m^, and some hypergeometric functions had a region of convergence which allowed us to 
study such integrals in that point. On the other hand, such hypergeometric series had a direct analytic continuation 
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into another one which did not allow us to consider s = 4m^. The answer is in when we carry out an analytic 
continuation and in the process cross a branch cut, this analytic continuation is not unique and poles do appear. In 
that case, there were simple and double poles, and we had two possible cases: F3 Y,F2 and F3 H2. In our present 
case, {MV)ii and (7V1'P)9 have third order poles, and [MV)'; is finite. It depends on the kinematical region we wish 
to study the pertinent integral. 

Back in the 50's and 60's |^ a lot of research was done in order to study singularities of Feynman integrals. One of 
the results we borrow from them is that a graph like the one considered in the present section can not have singularities 
in the physical sheet. This is a well-known theorem due to Eden (1952). So, in a full calculation of photon-photon 
scattering the poles must cancel out. 

In order to extract specific pole structures of these integrals, we can proceed just like in our previous works on the 
subject [ p^ . Expand gamma functions around e = and use Taylor expansion in the hypergeometric series, so that 
we are left with parametric derivatives of hypergeometric functions. The reader can see detailed calculations in the 
above referred paper; also in the appendix of pl| and in the next section. 



C. Non-Planar Double-Box integral 



Now we turn to massless non-planar double-box with six propagators, namely. 



(32) 



which represents the diagram of Figure 4. This diagram was also studied by Smirnov and Veretin that presented 
an explicit result in the case where all exponents were equal to minus one. On the other hand, we will write down 
results for arbitrary values of them. 
The generating functional 



d^qd^r exp{-ag2 - P{q - Pif -l{q + rf -0{q + r+p2f - (l>r^ - ij{r ~ p^f] 



2\D/2 



— ) exp 



-- {(3^ujs + a9u)t + (3e4>u) 



(33) 



can be integrated out without difficulty. We define A3 = a-|-/3 + 7 + and Q = a'^ + a9 + (3^ + (39 + \z{(t) + oj) ■ It is 
easy to see that Taylor expanding the above exponential will give us three series, while multinomial expansion for Q 
other twelve series. The equations that form the system to be solved come from the propagators, six altogether, and 
an additional constraint originates from the multinomial expansion. 

So, at the end of the day we are left with (15 — 7) = 8-fold series. Their variables are either. 



t s t 

~ : T ' ' 
S I U 



u u s 
t s u 



and/or unity. The simplest hypergeometric series representations for AfVe are double series, 

(-*|X2)(-m|X3)(-a'|X23)(-Z|X23) i-t/s)^- {u/s)^^ 



(1 + fc - C7'\X23){1 + n- a'\X3){l+j - a'\X2) X2lX3\ 



(34) 



(35) 



where a' — i+j + k + l + m + ji + D, is the sum of exponents and dimension and 

= {^j\<j'){-k\a'){-n\cr'){a' + D/2\ - 2a' - D/2){i + j+m + n + D\-i-j- D/2) 
X {k + l + m + n + D\-m-n~ D/2){i + i + k + l + D\~k~l~ D/2). 

The above hypergeometric series reduces to an Appel F2 function in the special case where k — —1^ 



—t u 
s s 



(36) 



(37) 



If we take all the exponents to be equal to minus one, we have 
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_gr(6 ^ D)V^{D/2 - 2)r3(D - 5) 
T(3lV2~~6)r^(I?~"4) 



E 



(1|X2)(1|X3)(1|X23) 



(6-i?|X3)(6-i?|X2 



D z,_6 r(6-i^)r^p/2-2)r^(i^-5) 

r(3i:>/2-6)r3(i:)-4) ^ 



1,1,1 

6 -D, 6-13 



s s 



(38) 
(39) 

This series is the Appel |15[] F2 hypergeometric function which converges when \t/ s\ < 1, \u/s\ < 1 and \t/s\ + \u/s\< 
1. Note that the above resuh has a double pole, 1/e^, just like in Tausk's Q and in Smirnov and Veretin's [Q. 
works. 

One could now ask: Would the divergent pieces come from a few terms in the series or from all of them? Clearly, 
in the present case, divergent factors which generate double and simple poles come from the gamma functions. To 
write down these poles explictly we have to Taylor expand also the hypergeometric function F2 , then we have around 
e = (£» = 4-2e). 



' 2,2 



r(3i:»/2-6)r3(i:>-4) 



1,1,1 
6-D,6-D 



t u 
s s 



(40) 



t u 

S S 



where 7^ is the Euler's constant and 



d^F2 = —F2 





-'-'-) 


V 7,7' 


s s J 



a=/3=/3'=l, 7=7'=2 



are called parametric derivatives of hypergeometric functions and can be calculated using Davydychev's pq ] 
approach. (r is given below.) 

In order to rewrite the above result as polilogarithm functions, such as Li2, Lis and usual logarithms we have to 
use an integral representation for i^2, 



r(/3)r(/30r(7 - /3)r(y - p') 
r(7)r(7') 



Fo 



7,7' 



— / dxidx2 



(1 - XiPi - X2P2Y 



(41) 



where Xi + X2 — 1. This task in not an easy one at all, since the second derivatives give a very cumbersome result 
and direct comparison between our result and Tausk's |^ is not possible (we were not able to show that both results 
are equivalent analytically). 

However, the important special case of forward scattering can not be read off directly from Tausk's result. On the 
contrary, making i = or u = is extremely simple in our eq.(|3^). Let m = 0, then one of the series (in X3 index) 
does not contribute and 



1,1,1 
6-i:i,6-D 



-,o 

s 



2F1 



1,1 
D 



r(6 - L')r(4 - £>) 



(42) 



since s + t + u — Q we were able to sum up the series 2^"!, so d jd^ can be rewritten in terms of d jdt and parametric 
derivatives are easily calculated from eq.(E2h. We quote the result up to order e". 



D „-2-2e 



201 204 

2? " 



83 + leV- (1) - 16V' (2) 



+4 



-134 



35?A"(1) -35V'"(2) +0(e) | 



(43) 



where '0(z)r(z) — r'(z), ip'{z), ij/'iz) their derivatives and 

^_ r{5-D)T^{D/2-l) 
T{3D/2-5){5- D) ' 

which in fact can be calculated to arbitrary order, since we start from an exact result. 
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Observe that the result has poles of higher orders, namely and e~^, which come from the fact we have taken 
u = 0, which also mean that the result has a branch point in u = 0, a well-known fact from the theory of hypergeometric 
functions. 

Of course, the 7x15 system of linear algebraic equations defined by this integral, the 7x7 solvable subsystems 
are divided into the following categories: among the grand total of 6,435 possible solutions, 3,519 have no solution 
at all. Among the remaining 2, 916 relevant solutions, NDIM provides other kind of series, such as 7-fold series and 
5-fold series. And all of them have symmetries among s,t and u, namely, 

(P3 P4j j ^ i m, t u) , {p2 ^ P3, I ^ n, k ^ m, t s) , (p2 ^ P4, j I, i ^ k, s ^ u) , (44) 

So, for each hypergeometric series representations provided by NDIM there are other two, also originated from the 
system of algebraic equations, which represent the same integral and can be transformed into the first using (|44|). 
This is the case of (p5[), and it is stated also by Tausk ||], i.e., the diagram is completely symmetric (through analytic 
continuation) under external legs permutations. 




FIG. 3. Scalar massless non-planar double-box with six propagators. The labels in the internal lines represent the exponents 
of propagators, see (^. 

Despite the complicated form of non-planar double-box with six propagators, the result we obtained is very simple, 
a double hypergeometric series, which in the special case of fc = — 1 reduces to F2 Appel function. Tausk presented a 
result for the same graph in terms of di- and trilogarithms, in the special case where all exponents are equal to minus 
one. If we recall [Q, Davydychev calculated a scalar integral for photon-photon scattering and transformed the four 
F2 Appel hypergeometric functions into dilogarithnis. The same occurs here, in one hand di and trilogarithms and in 
the other double hypergeometric series, namely F2. 



III. CONCLUSION 



In this work we considered covariant gauge scalar pentabox and non-planar double-box integrals. In the former we 
considered two cases: all virtual particles being massless, and in the other, six of them having the same mass while 
the seventh is massless. The latter was calculated in the massless case and arbitrary exponents of propagators, a 
result which was missing in the literature. No reduction formulas, i.e., rules to connect a given diagram with simpler 
ones, were used. The results are given in terms of hypergeometric series and in different kinematical regions. 
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